IFUP-TH/2000-02 

The critical equation of state of three-dimensional XY systems 



Massimo Campostrini ^, Andrea Pelissetto^, Paolo Rossi ^, Ettore Vicari^ 
^ Dipartimento di Fisica delVUniversitd di Pisa and I.N.F.N., 1-56126 Pisa, Italy 
^ Dipartimento di Fisica delVUniversitd di Roma I and I.N.F.N., 1-00185 Roma, Raly 
e-mail: campoOmailbox . dif i .unipi . it, rossiOmailbox . dif i .unipi . it, 
Andrea. PelissettoOromal . inf n. it, vicariOmailbox . dif i .unipi . it 

(February 1, 2008) 

Abstract 

We address the problem of determining the critical equation of state of 
three-dimensional XY systems. For this purpose we first consider the small- 
field expansion of the effective potential (Helmholtz free energy) in the high- 
temperature phase. We compute the first few nontrivial zero-momentum n- 
point renormalized couplings, which parametrize such expansion, by analyz- 
ing the high-temperature expansion of an improved lattice Hamiltonian with 
suppressed leading scaling corrections. 

These results are then used to construct parametric representations of the 
critical equation of state which are valid in the whole critical regime, satisfy 
the correct analytic properties (Griffith's analyticity) , and take into account 
the Goldstone singularities at the coexistence curve. A systematic approxi- 
mation scheme is introduced, which is limited essentially by the number of 
known terms in the small-field expansion of the effective potential. 

From our approximate representations of the equation of state, we derive 
estimates of universal ratios of amplitudes. For the specific-heat amplitude ra- 
tio we obtain /A~ = 1.055(3), to be compared with the best experimental 
estimate A+ /A" = 1.054(1). 

Keywords: Critical Phenomena, three-dimensional XY systems, A- 
transition of ^He, Critical equation of state. Effective potential, High- 
Temperature Expansion. 
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I. INTRODUCTION 



In the theory of critical phenomena, continuous phase transitions can be classified into 
universality classes determined only by a few basic properties characterizing the system, such 
as the space dimensionality, the range of interaction, the number of components and the 
symmetry of the order parameter. The renormalization-group theory predicts that, within 
a given universality class, the critical exponents and the scaling functions are the same for 
all systems. Here we consider the XY universality class, which is characterized by a two- 
component order parameter and effective short-range interactions. The lattice spin model 
described by the Hamiltonian 

Hl = - J ^ Si - Sj + '^hi- Si, (1) 

<ij> i 

where Si is a two-component spin satisfying Sj ■ Sj = 1, is one of the systems belonging to the 
XY universality class. It may be viewed as a magnetic system with easy-plane anisotropy, 
in which the magnetization plays the role of order parameter and the spins are coupled to 
an external magnetic field h. 

The superfiuid transition of ^He, occurring along the A-line Tx{P) (where P is the pres- 
sure), belongs to the three-dimensional XY universality class. Its order parameter is related 
to the complex quantum amplitude of helium atoms. It provides an exceptional opportunity 
for an experimental test of the renormalization-group predictions, thanks to the weakness of 
the singularity in the compressibility of the fiuid and to the purity of the sample. Moreover, 
experiments in a microgravity environment lead to a reduction of the gravity-induced broad- 
ening of the transition. Recently a Space Shuttle experiment performed a very precise 
measurement of the heat capacity of liquid helium to within 2 nK from the A-transition 
obtaining an extremely accurate estimate of the exponent a and of the ratio A^/A^ of the 
specific-heat amplitudes: 

a = -0.01285(38), = 1.054(1). (2) 

These results represent a challenge for theorists because the accuracy of the test of the 
renormalization-group prediction is now limited by the precision of the theoretical calcula- 
tions. We mention the best available theoretical estimates for a: a = —0.0150(17) obtained 
using high-temperature expansion techniques 0, a = —0.0169(33) from Monte Carlo simu- 
lations using finite-size scaling techniques 0, a = —0.011(4) from field theory [Q. The close 
agreement with the experimental data clearly supports the standard renormalization-group 
description of the A-transition [Q. 

In this paper we address the problem of determining the critical equation of state char- 
acterizing the XY universality class. The critical equation of state relates the thermody- 
namical quantities in the neighborhood of the critical temperature, in both phases. It is 
usually written in the form (see e.g. Ref. |^) 

H = MM^-^f{x), xoctM-^/^, (3) 

where f{x) is a universal scaling function (normalized in such a way that /(—I) = and 
/(O) = 1). The universal ratios of amplitudes involving quantities defined at zero- momentum 
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(i.e. integrated in the volume), such as specific heat, magnetic susceptibihty, etc., can be 
obtained from the scaling function f{x). 

It should be noted that, for the A-transition in "^He, Eq. is not directly related to the 
conventional equation of state that relates temperature and pressure. Moreover, in this case, 
the field H does not correspond to an experimentally accessible external field, so that the 
function appearing in Eq. (|^) cannot be determined directly in experiments. The physically 
interesting quantities are universal amplitude ratios of quantities formally defined at zero 
external field. 

As our starting point for the determination of the critical equation of state, we com- 
pute the first few nontrivial coefficients of the small-field expansion of the effective poten- 
tial (Helmholtz free energy) in the high-temperature phase. For this purpose, we analyze 
the high-temperature expansion of an improved lattice Hamiltonian with suppressed lead- 
ing scaling corrections 0,0. If the leading non-analytic scaling corrections are no longer 



present, one expects a faster convergence, and therefore an improved high-temperature ex- 
pansion (IHT) whose analysis leads to more precise and reliable estimates. We consider a 
simple cubic lattice and the (p^ Hamiltonian 

n= -pJ2$^-^y+E[fl + H$l-'^y]^ (4) 

{x,y) ^ 

where (x, y) labels a lattice link, and (p^ is a real two-component vector defined on lattice 
sites. The value of A at which the leading corrections vanish has been determined by Monte 
Carlo simulations using finite-size techniques obtaining A* = 2.10(6). In Ref. [Q] we 
have already considered the high-temperature expansion (to 20th order) of the improved 
0^ Hamiltonian (^ for the determination of the critical exponents, achieving a substantial 
improvement with respect to previous theoretical estimates. IHT expansions have also been 



considered for Ising-like systems [|I0|, obtaining accurate determinations of the critical ex- 
ponents, of the small-field expansion of the effective potential, and of the small-momentum 
behavior of the two-point function. 

We use the small-field expansion of the effective potential in the high-temperature phase 
to determine approximate representations of the equation of state that are valid in the 
whole critical regime. To reach the coexistence curve {t < 0) from the high-temperature 



phase {t > 0), an analytic continuation in the complex t-plane [£|,11| is required. For this 
purpose we use parametric representations |T2|-|T^], which implement in a rather simple way 
the known analytic properties of the equation of state (Griffith's analyticity). This approach 
was successfully applied to the Ising model, for which one can construct a systematic ap- 



proximation scheme based on polynomial parametric representations [|11|] and on a global 
stationarity condition [0. This leads to an accurate determination of the critical equation 
of state and of the universal ratios of amplitudes that can be extracted from it |Tl],^,p!0|. XY 



systems, in which the phase transition is related to the breaking of the continuous symme- 
try 0(2), present a new important feature with respect to Ising-like systems: the Goldstone 
singularities at the coexistence curve. General arguments predict that at the coexistence 
curve (t < and H 0) the transverse and longitudinal magnetic susceptibilities behave 
respectively as 

Xr=^. x, = — ^H''^-\ (5) 
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In our analysis we will consider polynomial parametric representations that have the correct 
singular behavior at the coexistence curve. 

The most important result of the present paper, from the point of view of comparison 
with experiments, is the specific-heat amplitude ratio; our final estimate 

A^jA- = 1.055(3) (6) 

is perfectly consistent with the experimental estimate (^, although not as precise. 

The paper is organized as follows. In Sec. |I| we study the small-field expansion of the 
effective potential (Helmholtz free energy). We present estimates of the first few nontriv- 
ial coefficients of such expansion obtained by analyzing the IHT series. The results are 
then compared with other theoretical estimates. In Sec. |T|, using as input parameters the 
critical exponents and the known coefficients of the small-field expansion of the effective 
potential, we construct approximate representations of the critical equation of state. We 
obtain new estimates for many universal amplitude ratios. These results are then compared 
with experimental and other theoretical estimates. 



II. THE EFFECTIVE POTENTIAL IN THE HIGH-TEMPERATURE PHASE 



A. Small-field expansion of the effective potential in the high-temperature phase 

The effective potential (Helmholtz free energy) is related to the (Gibbs) free energy of 
the model. If M = (0) is the magnetization and H the magnetic field, one defines 

J^{M) = M-H-^\ogZ{H), (7) 

where Z{H) is the partition function and the dependence on the temperature is always 
understood in the notation. 

In the high-temperature phase the effective potential admits an expansion around M = 0: 

oo 1 

= HM) - ^(0) = J: J^_^2jM''- (8) 
This expansion can be rewritten in terms of a renormalized magnetization ip 

A^=imV + E^'"^(^^2y^' (9) 

where 

t is the reduced temperature, x and ^ are respectively the magnetic susceptibility and the 
second-moment correlation length 
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x = Y.^Kmc.{x)), (11) 

X 

and m = 1/^. In field theory ip is tlie expectation value of the zero-momentum renormalized 
field. The zero-momentum 2j-point renormalized constants g2j approach universal constants 
(which we indicate with the same symbol) for t — 0. By performing a further rescaling 

^ = (12) 

in Eq. (|^), the free energy can be written as 

777 

^T=—A{z), (13) 
9i 



where 



Aw4.= + i.' + Ep^r,,.=^ (14) 



and 



-2, = J > 3. (15) 

One can show that z oc t~^M, and that the equation of state can be written in the form 

H oc (16) 

(7Z 



B. Zero-momentum renormalized couplings by IHT expansion 

To compute the high-temperature series of the four-point coupling and of the effective- 
potential parameters r2j, we rewrite them in terms of the zero-momentum connected 2j-point 
Green's functions X2j 

X2j = {(f>ai{0)(f)aAx2)--(j)a,{x2j-l)(f)a,{x2j))c (17) 

X2,---,X2j 

ix = X2)- For generic A^- vector models we have 



* = -—2lse '''' 



and 
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3(iV + 4) xl ' ^ ^ 



280(iV + 2)x6X2 , 35(iV + 2)2 xsxl 



^8 = 280 - —- -r^-^ + 



3(iV + 4) xl 9(iV + 4)(iV + 6) xl 

7700(iV + 2) X6X2 ^ 350(iV + 2)^ ^^^^ 
rio = 15400 5 — h 



(iV + 4) xl (iV + 4)2 xl 
1400(iV + 2)2 xsxi 35(iV + 2)3 xioX^ 



3(iV + 4)(Ar + 6) xl 3(Ar + 4)(Ar + 6)(iV + 8) xl ' 

The formulae relevant for the XY universality class are obtained setting N = 2. 

Using the 0*^ lattice Hamiltonian (^), we have calculated x and m2 = J2x^'^{(p{0)'P{^)) 
to 20th order, X4 to 18th order, Xe to 17th order, xs to 16th order, and Xio to 15th order, 
for generic values of A. The IHT expansion, i.e. with suppressed leading scaling corrections, 
is achieved for A = 2.10(6) [If. In Table | we report the series of m2, X; X4; Xe, Xs and 
Xio for A = 2.10. Using Eqs. (|18D and (pJ]) one can obtain the HT series necessary for the 
determination of and r2j. We analyzed the series using the same procedure applied to 
the improved high-temperature expansions of Ising-like systems in Ref. |iy]. In order to 



estimate the fixed-point value of and r2j, we considered Fade, Dlog-Pade and first-order 
integral approximants of the series in (3 for A = 2.10, and evaluated them at (3c- We refer to 
Ref. [|10| for the details of the analysis. Our estimates are 

^4 = 21.05(3 + 3), (20) 

^ = ^^74 = 1.396(2 + 2), (21) 

re = 1.951(11 + 3), (22) 
rg = 1.36(6 + 3). (23) 

We quote two errors: the first one is related to the spread of the approximants, while 
the second one gives the variation of the estimate when A varies between 2.04 and 2.16. 
In addition, we obtained a rough estimate of rio, i.e. rio = —13(7). For comparison, we 
anticipate that the analysis of the critical equation of state using approximate parametric 
representations will lead to the estimate rio = —10(3). From the estimates of gl and 
one can obtain corresponding estimates for the zero-momentum renormalized couplings, 
g2j = r2jgi~^ with j > 2. 

Table |T| compares our results (denoted by IHT) with the estimates obtained using other 
approaches, such as the high-temperature expansion (HT) of the standard lattice spin model 
(HI) fTsHlTl, field-theoretic methods based on the fixed- dimension d = 3 gf-expansion 0,|T8|j2O 



and on the e-expansion p6|j21|j2^ . The fixed-dimension field-theoretic estimates of (74 have 
been obtained from the zero of the Callan-Symanzik /3-function, whose expansion is known 
to six loops |T9|. In the same framework g^ and gs have been estimated from the analysis 



of the corresponding four- and three- loop series respectively |18]. The authors of Ref. | T8[ 



argue that the uncertainty on their estimate of ge is approximately 0.3%, while they consider 
their value for g^ much less accurate. The e-expansion estimates have been obtained from 
constrained analyses of the four-loop series of (74 and the three- loop series of r2j. 
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TABLE I. Coefficients of the high-temperature expansion of m2, X4) XGi XSj Xio- They have 
been obtained using the Hamiltonian (^) with A = 2.10. 



t 


ni2 


X 


X4 








0.82305062235187783163838 


-0.18361621025768068492172 


1 


1.01611849043072013077966 


2.03223698086144026155932 


-1.81350523351772810305606 


2 


5.01790173559352893587925 


4.41339999108761956819390 


-10.2403328491873632045378 


3 


17.1372887645210686867092 


9.49465465075720590174531 


-44.8928073800743626690682 


4 


50.5729609164717064337798 


19.8432084757850795073880 


-168.753294310535712758532 


5 


136.975729153185634797939 


41.3105316077169269675024 


-572.899325351072412028155 


6 


351.547955204015543873533 


84.8657665346587543517940 


-1806.05498365521624417337 


7 


867.882609960590948513531 


173.919877080583932668034 


-5384.52486767607368575700 


8 


2082.72988964511892008013 


353.805866554914925628451 


-15362.4317708302607812112 


9 


4887.89348922134368695104 


718.544573182206654200415 


-42306.4841196583464482445 


10 


11271.6038365544561486584 


1452.55696588171500306095 


-113150.461045583807761395 


11 


25618.1126339930289598217 


2932.78243299692312158659 


-295305.926162214852952064 


12 


57532.0447298988841 107389 


5902.70244465199387875428 


-754781.109517718395587745 


13 


127889.749773305781861722 


11869.0080781568501056056 


-1894806.32404641680977111 


1 A 
14 




zooil.UDzzloUUizUbozoUzz ( 


— 4ooz /oy.D44ol / /zoozo/oy41 


15 


616366.979493029224503143 


47733.2907485305100325248 


— 11414677.5171522019349233 


16 


1339133.76964227829126980 


95522.7109197203117005540 


—27486804.4332023911210055 


17 


2892420.77412897864771425 


191043.098791148001757550 


—65472437.4333935066305991 


1 O 

18 


6215040.60947250914305650 


OOlCT/^r* i^C\ A C\'\ C\i^ A C\r' A '\ '1C\ O O y1 /I 

381560.204019249641798944 


— 154435238.761288809747194 


19 


13292241.3642748890173429 


761691.955581673385411208 




20 


28309542.5693146536002667 


1518865.25572011244933658 




i 


X6 


X8 


Xio 





0.41651792496748423277383 


-2.07276872395008612032057 


17.8224779282297038884825 


1 


8.19357083273717982769144 


-66.6409389277329689706620 


845.247991215093077905035 


2 


81.6954441314190452820070 


-1011.50102952099242860043 


17983.4651630934745766842 


3 


577.734561420010600288632 


-10327.4974814820902607880 


248053.171125799025029980 


4 


3288.50395189429008814810 


-81557.4477403308614541401 


2572480.03564061669203194 


5 


16096.6021109696857145564 


-536525.345841149579660234 


21718902.5365937343351130 


6 


70442.2092664445625747715 


-3075097.70592321924332973 


156737903.631432308532083 


7 


282678.282078376456985396 


-15817337.6236306593687791 


998806344.593851737549288 


8 


1058442.95252211500336034 


-74543191.4965166803297285 


5750852421.81974869997198 


9 


3744677.84665924694446919 


-326778866.937610399640322 


30428100948.8992362353297 


iU 


izDOO /OU.o /0ozy4zo r oUio 


— io4rozoior .yyor /o4odUzoo4 


i4yoD00iDZ0Z.i4U44i0U44 ( 4 


11 


40959287.0765118834279901 


-5277067971.27466652475816 


694034037396.324618345020 


12 


128268696.512235415964465 


-19750695382.3283907899788 


3046464352917.44305602911 


13 


389824946.322878356697650 


-71066114622.3163753581174 


12757807364908.6407194560 


14 


1153973293.83499223649654 


-246972985495.952557381976 


51244942192011.6695801488 


15 


3337461787.55665279189948 


-832179913942.254869433972 


198319718730601.629272462 


16 


9454326701.58705373223814 


-2727546459564.59095483523 




17 


26288539415.6469228911099 







7 



TABLE II. Estimates of g = 554/(247r), re, and rg. 





IHT 


HT 


d = 3 g-exp. 


e-exp. 


g 


1.396(4) 


1.411(8), 1.406(8) (1|1 


1.403(3) i 


1.425(24) l22,16| 






1.415(11) PI 


1.40 m 




re 


1.951(14) 


2.2(6) 


1.967 |T8| 


1.969(12) l2^j21| 




1.36(9) 




1.641 |18| 


2.1(9) |2|,21| 



III. THE CRITICAL EQUATION OF STATE 

A. Analytic properties of the scaling equation of state 

From the analysis of the IHT series we have obtained the first few non-trivial terms of the 
small-field expansion of the effective potential in the high-temperature phase. This provides 
corresponding information for the equation of state 

Hoct'^^F{z), (24) 

where z oc Mt~^ and, using Eq. (IT 



F{z) = ^ = ^ + 1^' + E F2m-iz'"^~' (25) 

m=3 



with 



2m- 1 



(2m - 1)! 



r2m- (26) 



The function H{M^t) representing the external field in the critical equation of state (|24D 
satisfies Griffith's analyticity: it is regular at M = for t > fixed and at t = for M > 
fixed. The first region corresponds to small z in Eq. (^4]), while the second is related to large 
z, where F{z) has an expansion of the form 

F{z) = z^Y. F^z-'''^. (27) 

n=0 

To reach the coexistence curve, i.e. t < and H = 0, one should perform an analytic 
continuation in the complex t-plane []8|,pA| . The spontaneous magnetization is related to the 



complex zero Zq of F{z). Therefore, the description of the coexistence curve is related to 
the behavior of F{z) in the neighbourhood of 2:0- 



B. Goldstone singularities at the coexistence curve 

The physics of the broken phase of A^-vector models (including XY systems which cor- 
respond to = 2) is very different from that of the Ising model, because of the presence 
of Goldstone modes at the coexistence curve. The singularity of xl for t < and H ^ 
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is governed by the zero-temperature infrared- stable fixed point p3|-p5|. This leads to the 
prediction 



fix) ^Cf{l + 



for 



-1, 



(28) 



where x oc tM~^^^ and f{x) is the scaling function introduced in Eq. (0) (as usual, x = —1 
corresponds to the coexistence curve). This behavior at the coexistence curve has been 
verified in the framework of the large- iV expansion to 0{1/N) (i.e. next-to-leading order) 
|2|,|2|. 

The nature of the corrections to the behaviour (^) is less clear. Setting u = 1 + x and 
y = HM~^, it has been conjectured that u has the form of a double expansion in powers of 
y and 7/('^~2)/2 j^gg^j^ |^]-^g coexistence curve ||27|j28|j25[| , i.e. for y 







uj = l + x = ciy + C2y^~'^^ + diy' + d2y'''''' + d^y'~' + 



2-e/2 



(29) 



where e 
29| 



SIS 

,1/2 



= A — d. This expansion has been derived essentially from an e-expansion analy- 
Note that in three dimensions this conjecture predicts an expansion in powers of 
y"'", or equivalently an expansion of f{x) in powers of uj for a; — > 0. 

The asymptotic expansion of the d-dimensional equation of state at the coexistence curve 



has been computed analytically in the framework of the large- iV expansion |22[, using the 
0{1/N) formulae reported in Ref. It turns out that the expansion (^Uj) does not strictly 
hold for values of the dimension d such that 



2 < = 2 



2m 



<4, 



n 



for < m < n, m, G IN. 



In particular, in three dimensions one finds |^ 



UJ 



1 



l + j-{f,{uj)+\0gLuf2iuj)) + 0iN-'] 



(30) 



(31) 



where the functions fi{uj) and /2(t^) have a regular expansion in powers of u. Moreover, 

Uu;)=0{uj'), (32) 

so that the logarithms affect the power expansion only at the next-next-to-leading order. A 
possible interpretation of the large- analysis is that the expansion ( pT)) holds for all values 
of A^, so that Eq. (E^ is not correct due to the presence of logarithms. The reason of their 
appearance is however unclear. Neverthless, it does not contradict the conjecture that the 
behavior near the coexistence curve is controlled by the zero-temperature infrared-stable 
Gaussian fixed point. In this case logarithms would not be unexpected, as they usually 
appear in the reduced-temperature asymptotic expansion around Gaussian fixed points (see 
e.g. Ref. [0). 



C. Parametric representations 

In order to obtain a representation of the critical equation of state that is valid in the 
whole critical region, one may use parametric representations, which implement in a simple 

9 



way all scaling and analytic properties |T^-|T4[. One may parametrize M and t in terms of 
R and 9 according to 

M = moR^m{e), (33) 

t = R{l-0^), 
H = hoRf^^hie), 

where ho and mg are normalization constants. The variable R is nonnegative and measures 
the distance from the critical point in the (t, H) plane; it carries the power-law critical 
singularities. The variable 6 parametrizes the displacements along the line of constant R. 
The functions m(6') and h{9) are odd and regular at 6* = and at 6^ = 1. The constants tuq 
and ho can be chosen so that m{6) = 6 + 0{6^) and h{6) = 6 + 0{d^). The smallest positive 
zero of h{6), which should satisfy 6*0 > 1, represents the coexistence curve, i.e. T < T^. and 
H ^0. 

The parametric representation satisfies the requirements of regularity of the equation 
of state. Singularities can appear only at the coexistence curve (due for example to the 
logarithms discussed in Sec. |111BD , i.e. for 9 = 6q. Notice that the mapping ( PHj ) is not 



invertible when its Jacobian vanishes, which occurs when 



Y{e) = {1- e^)m'{e) + 2i3em{e) = o. (34) 

Thus the parametric representations based on the mapping (^) are acceptable only if < 6i 
where 6i is the smallest positive zero of the function Y{6). One may easily verify that the 
asymptotic behavior (^) is reproduced simply by requiring that 

h{e) ^ - ef for e e^. (35) 

The relation among the functions 171(6), h{9) and F{z) is given by 

z = pm{e){l-e^y\ (36) 

Fizie)) = p(i-e'y''hie), (37) 



where p is a free parameter |TT],|TO[. Indeed, in the exact parametric equation the value of p 
may be chosen arbitrarily but, as we shall see, when adopting an approximation procedure 
the dependence on p is not eliminated. In our approximation scheme we will fix p to ensure 
the presence of the Goldstone singularities at the coexistence curve, i.e. the asymptotic 
behavior (|35|). Since z = p6 + 0{9^), expanding m{9) and h{6) in (odd) powers of 9, 



m 



(^)=^ + E"^2„+i^'"+\ (38) 



n=l 



n=l 



and using Eqs. (pq ) and (pTl), one can find the relations among p, m2n+i, ^2n+i and the 
coefficients F2„+i of the expansion of F{z). 

One may also write the scaling function f{x) in terms of the parametric functions m{6) 
and h{ey. 
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0'^ 

^0 



fix) 



m{e) 
m(l) 



m{9) 
1 -s 



m 



In App. we report the definitions of some universal ratios of amphtudes that have been 
introduced in the hterature, and the corresponding expressions in terms of the functions m{9) 
and h{e). 



D. Approximate polynomial representations 

In order to construct approximate parametric representations we consider polynomial 
approximations of m{6) and h{6). This kind of approximation turned out to be effective in 
the case of Ising-like sistems [p]T|p!0| . The major difference with respect to the Ising case is 
the presence of the Goldstone singularities at the coexistence curve. In order to take them 
into account, at least in a simplified form which neglects the logarithms found in Eq. (Plj), 
we require the function h{6) to have a double zero at as in Eq. (pSf). Polynomial schemes 
may in principle reconstruct also the logarithms, but of course only in the limit of an infinite 
number of terms. 

In order to check the accuracy of the results, it is useful to introduce two distinct schemes 
of approximation. In the first one, which we denote as (A), h{6) is a polynomial of fifth 
order with a double zero at and m{6) a polynomial of order (1 + 2n): 

scheme (A) : m{e) =o(^l + J2 c*^^*) > (40) 

h{e) = e{i-eyeiy. 

In the second scheme, denoted by (B), we set 

scheme (B) : m{9) = 9, (41) 

hi9) =9{l- 9'/9iy (^1 + Ec.^'') • (42) 

Here h{9) is a polynomial of order 5 + 2n with a double zero at ^o- In both schemes the 
parameter p is fixed by the requirement (|55D, while 9o and the n coefficients q are determined 
by matching the small-field expansion of F{z). This means that, for both schemes, in order 
to fix the n coefficients q we need to know n + 1 values of r2j, i.e. rg, ..■rQ+2n- Note that for 
the scheme (B) 

Y{9) = 1-9"^ + 2(39^ (43) 

independently of n, so that = (1 — 2f3)~^. Concerning the scheme (A), we note that the 
analyticity of the thermodynamic quantities for \9\ < 9o requires the polynomial function 
Y{9) not to have complex zeroes closer to origin than 9o. 

In App. ^ we present a more general discussion on the parametric representations. 
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TABLE III. Results for the parameters and the universal amplitude ratios using the scheme 
(A), cf. Eq. (|0D, and the scheme (B), cf. Eq. @). The label ab ove each column indicates the 
scheme, the number of terms in the corresponding polynomial, and the input parameters employed, 
beside the critical exponents a and rj. Note that the quantities reported in the first four lines do 
not have a physical meaning, but are related to the particular parametric representation employed. 
Numbers marked with an asterisk are inputs, not predictions. 





[{J\) n — i, rg, rgj 


[(B) n = 1; re, rgj 


[(B) n — 2; rg, rg, rioj 


\(T3\ „ 0- ly rn 4+ /Zl^l 

[{H} n — Z, rg, rg, A /A \ 


p 


2.22(3) 


2.07(2) 


2.04(5) 


2.01(4) 


ol 




3.84(10) 


2.97(10) 


2.8(4) 


2.5(2) 


Cl 


-0.024(7) 


0.28(3) 


0.10(6) 


0.15(5) 


C2 








0.01(2) 


0.02(1) 


no 


-9.6(1.1) 


-11(2) 


*-13(7) 


-7(5) 


A+/A- 


1.055(3) 


1.057(2) 


1.055(3) 


*1.054(1) 


Rc 


0.123(8) 


0.113(3) 


0.118(7) 


0.123(6) 


R4 


7.5(3) 


7.9(2) 


7.8(3) 


7.6(2) 


R* 


0.353(3) 


0.350(2) 


0.352(3) 


0.354(3) 




1.38(9) 


1.51(3) 


1.47(8) 


1.41(6) 


7700 
^0 


0.0303(3) 


0.0301(3) 


0.0302(3) 


0.0304(3) 




5(4) 


62(41) 




15(10) 



E. Results 



As input parameters for the determination of the parametric representations, we use 
the best available estimates of the critical exponents, which are a = —0.01285(38) (from 
the experiment of Ref. |]l[), 1] = 0.0381(3) (from the high-temperature analysis of Ref. [0). 
Moreover we use the following estimates of r2j: rg = 1.96(2) which is compatible with all the 
estimates of rg reported in Table ||, and rg = 1.40(15) which takes somehow into account 
the differences among the various estimates. 

The case n = of the two schemes (A) and (B) is the same, and requires the knowledge 
of a, rj and rg. Unfortunately this parametrization does not satisfy the consistency condition 
9q < Of = (1 — 2(3)^^. Both schemes give acceptable approximations for n = 1, using rg 
as an additional input parameter. The numerical values of the relevant parameters and the 
resulting estimates of universal amplitude ratios (see the appendix for their definition) are 



shown in Table |ITI[ The errors reported are related to the errors of the input parameters 
only. They do not take into account possible systematic errors due to the approximate 
procedure we are employing. We will return on this point later. 

In Figs. |l| and |^ we show respectively the scaling functions F{z) and f{x), as obtained 
from the approximate representations given by the schemes (A) and (B) for n = 1, using the 
input values a = —0.01285, r] = 0.0381, rg = 1.96 and rg = 1.4. The two approximations of 
F(z) are practically indinstinguishable in Fig. ^ This is also numerically confirmed by the 
estimates of the universal costant (reported in Table p^ , which is related to the large- 2; 
behavior of F{z): 

F{z) ^ F^z^ for z^oo. (44) 
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z 

FIG. 1. The scaling function F{z) vs. z. The convergence radius of the small-z expansion is 

1/2 

expected to be |zo| = — 2.8. 

This agreement is not trivial since the small- 2; expansion has a finite convergence radius 
given by \zq\ = i?^^ ~ 2.8. Therefore, the determination of F{z) on the whole positive 
real axis from its small- 2; expansion requires an analytic continuation, which turns out to 
be effectively performed by the approximate parametric representations we have considered. 
We recall that the large-z limit corresponds to the critical theory t = 0, so that positive 
real values of z describe the high-temperature phase up to t = 0. Instead, larger differences 
between the approximations given by the schemes (A) and (B) for n = 1 appear in the 
scaling function f{x), especially in the region a; < corresponding to t < (i.e. the region 
which is not described by real values of z). Note that the apparent differences for x > are 
essentially caused by the normalization of /(x), which is performed at the coexistence curve 
X = —1 and at the critical point a; = requiring /(—I) = and /(O) = 1. Although the 
large-x region corresponds to small z, the difference between the two approximate schemes 
does not decrease in the large-x limit due to their slightly different estimates of (see 
Table |IT|). Indeed, for large values of x, f{x) has an expansion of the form 

fix) =x'<Y. /r^"'"^ (45) 

ra=0 

with = R-\ 

We also considered the case n = 2, using the estimate rio = —13(7). In this case the 
scheme (A) was not particularly useful because it turned out to be very sensitive to rio, 
whose estimate has a relatively large error. Combining the consistency condition 6q < 61 
(which excludes values of rio ^ —10 when using the central values of a, r], rg and rg) with 
the IHT estimate of rio, we found a rather good result for /A~ , i.e. A'^/A' = 1.053(4). 
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fix) 




FIG. 2. The scaling function f{x) vs. x. 



On the other hand, the results for the other universal amplitude ratios considered, such as 
i?^, i?c, Rii etc., although consistent, turned out to be much more imprecise than those 
obtained for n = 1, for example = 1.2(3). This fact may be explained noting that, except 
for the small interval —10 ^ rio ^ —9 (this interval corresponds to the central values of the 
other input parameters), the function Y{6), cf. Eq. (|3^) , has zeroes in the the complex plane 
which are closer to the origin than 6q. Therefore, the parametric function g2{0) related to 
the magnetic susceptibility (see App. [A 2| ), and higher-order derivatives of the free-energy, 
have poles within the disk |^^| < Oq. On the other hand, in the case n = 1, 6q was closer to 
the origin than the zeroes of Y{6) for the whole range of values of the input parameters. 

For these reasons, for n = 2, we present results only for the scheme (B). Combining the 
consistency condition 9q < 9i, which restricts the acceptable values of rio (for example using 
the central estimates of the other input parameters it excludes values |rio| ^ 12), with the 
IHT estimate rio = ~13(7), we arrive at the results reported in Table ^ (third column of 
data). The reported value has been obtained using rio ~ —9. 

The coefficients q, reported in Table pTT|, turn out to be relatively small in both schemes, 
and decrease rapidly, supporting our choice of the approximation schemes. The results 
of the various approximate parametric representations are in reasonable agreement. Their 
comparison is useful to get an idea of the systematic error due to the approximation schemes. 
There is a very good agreement for A^/A~, which is the experimentally most important 
quantity. Our final estimate is 

A+/A- = 1.055(3). (46) 
We mention that approximately one half of the error is due to the uncertainty on the critical 
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TABLE IV. Estimates of universal amplitude ratios obtained in different approaches. The 
e-expansion estimates of Rc and R^, have been obtained by setting e = 1 in the O(e^) series 



calculated 


in Refs. 


34]. 




IHT PR 


HT, LT 




d=3 exp. 


e-exp. 


experiments 


A+/A- 


1.055(3) 


1.08 


35 


1 


1.056(4) m 


1.029(13) |7 


1.054(1) i 
1.058(4) 
1.067(3) |3|] 
1.088(7) |40|1 


"t 


0.353(3) 


0.361(4), 0.362(4) || 


0.3606(20) Ugl 


0.36 H 




Rc 


0.12(1) 








0.123(3) g 


0.106 






1.4(1) 










1.407 





exponent a, which unfortunately can be hardly improved by present theoretical means. 
The approximation schemes (A) and (B) with n = 1 provide independent results for rio, 
leading to the estimate rio = —10(3), which is agreement with the IHT result rio = —13(7). 
The determination of Cf, cf. Eq. (p8|) , turns out to be rather unstable, indicating that the 
approximate parametric representation we have constructed are still relatively inaccurate in 
the region very close to the coexistence curve. The constant cj is very sensitive to the values 
of the coefficients r2j. Improved estimates of r2j would be important especially for Cf. 

Finally, to further check our results, we applied again the scheme (B) for n = 2, replacing 
rio with the precise experimental estimate = 1.054(1) as input parameter. In practice 

we fix the coefficient C2 in such a way to obtain the experimental estimate of A'^/A^. The 
idea is to use the quantities known with the highest precision to determine, within our 
scheme of approximation, the equation of state and the corresponding universal amplitude 
ratios. The results are reported in the last column of Table |I|. 

From the results of Table |ITT| we arrive at the final estimates 

R+ = iA+y/''^f+ = 0.353(3), (47) 
Rc ^ = 0.12(1), (48) 
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Rx = ^77^ = 1-4(1), (49) 



^4 = -^ = koP = 7.6(4), (50) 

= lim z-^F(z) = 0.0303(3), (51) 
0<Cf< 20. (52) 

In Table |I^ we compare our results with the available estimates obtained from other 
theoretical approches and from experiments (for a review see e.g. Ref. P3|). Our precision for 



A'^ /A~ is comparable with the estimate reported in Ref. [^, obtained in the field-theoretic 
framework of the minimal renormalization without e-expansion, which is a perturbative 
expansion at fixed dimension d = 3. The agreement with the experimental result of Ref. |jl| 
is very good. 
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F. Conclusions 



Starting from the small-field expansion of the effective potential in the high-temperature 
phase, we have constructed approximate representations of the critical equation of state 
valid in the whole critical region. We have considered two approximation schemes based 
on polynomial representations that satisfy the general analytic properties of the equation 
of state (Griffith's analyticity) and take into account the Goldstone singularities at the 
coexistence curve. The coefficients of the truncated polynomials are determined by matching 
the small-field expansion in the high-temperature phase, which has been studied by lattice 
high-temperature techniques. The schemes considered can be systematically improved by 
increasing the order of the polynomials. However, such possibility is limited by the number of 
known coefficients r2j of the small-field expansion of the effective potential. We have shown 
that the knowledge of the first few already leads to satisfactory results, for instance for the 
specific-heat amplitude ratio. Through the approximation schemes we have presented in this 
paper, the determination of the equation of state may be improved by a better determination 
of the coefficients r2j, which may be achieved by extending the high-temperature expansion. 
We hope to return on this issue in the future. 

Finally we mention that the approximation schemes which we have proposed can be 
applied to other A^-vector models. Physically relevant values are N = 3 and = 4. 



The case = 3 describes the critical phenomena in isotropic ferromagnets The case 
= 4 is interesting for high-energy physics: it should describe the critical behavior of 
finite-temperature QCD with two flavours of quarks at the chiral- symmetry restoring phase 
transition [^]. 



APPENDIX A: UNIVERSAL RATIOS OF AMPLITUDES 

1. Notations 

Universal ratios of amplitudes characterize the critical behavior of thermodynamic quan- 
tities that do not depend on the normalizations of the external (e.g. magnetic) field, order 
parameter (e.g. magnetization) and temperature. Amplitude ratios of zero-momentum quan- 
tities can be derived from the critical equation of state. We consider several amplitudes 
derived from the singular behavior of the specific heat 

cj, = A±|tr", (Ai) 

the magnetic susceptibility in the high-temperature phase 

X = lcn-\ (A2) 
the zero-momentum four-point connected correlation function in the high temperature phase 

X4 = lctt-^-''\ (A3) 
the second-moment correlation length in the high-temperature phase 
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e = rt-\ (A4) 

and the spontaneous magnetization on the coexistence curve 

M = B\tf. (A5) 

Using the above normahzations for the amphtudes, the zero-momentum four-point couphng 
g4, cf. Eq. (|18|), can be written as 

In addition, one can also define amphtudes along the critical isotherm, such as 

XL = C'\H\-^^f^^. (A7) 



2. Universal ratios of amplitudes from the parametric representation 

In the following we report the expressions of the universal ratios of amplitudes in terms 
of the parametric representation ( ^3]) of the critical equation of state. 
The singular part of the free energy per unit volume can be written as 

J^ = homoR^~''9{e), (A8) 

where g{6) is the solution of the first-order differential equation 

(1 - e')g'ie) + 2(2 - a)egie) = Y{e)h{e) (A9) 



that is regular at 6 = 1. The function Y{9) has been defined in Eq. (|34D. The longitudinal 
magnetic susceptibility can be written as 

. . 2mm.il- erne) ^ 

mo Y{9) 

The function g2{0) must vanish at 9q in order to reproduce the predicted Goldstone singu- 
larities, according to 

g2{0) ^00-6 for 6^ 00. (All) 

From Eq. ( |A1CI|) we see that g2{0) satisfies this condition if h{9) ~ (^o — 0)^ for 6 Oq. 

From the equation of state one can derive universal amplitude ratios of quantities defined 
at zero momentum, i.e. integrated in the volume. We consider 

A^/A- = {91 - If-^y (A12) 

9{9o) 

^ = -a{l - a)(2 - aM - l?^[m{9o)]-'g{0), (A13) 

i?4 = = koP = [mi9o)r {91 - l)-'' , (A14) 

Rx = = (^^0 - ir'[m{9o)r'[m{l)]-'h{l), (A15) 
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Using Eqs. (pB]) and (^) one can compute F{z) and obtain the small-z expansion coef- 
ficients of the effective potential r2j. The constant F^, which is related to the behavior of 
F{z) for z ^ oo, cf. Eq. (0), is given by 

F^ = lim z-^F{z) = p^-^[m{l)]-^h{l). (A16) 

Using the relations (^) concerning the scaling function f{x), one can easily obtain the 
constant C/, which is related to the behavior of f{x) for x —1, cf. Eq. (^), 

Cf = lira fi + x)-^ f{x). (A17) 

We consider also the universal amplitude ratio = (A'^Y^^f^ which can be obtained from 
the estimates of i?4, Rc and (^4: 

4 ^ (A+yr^ /+ = (^^^ . (AI8) 

We mention that in the case of the superfluid helium it is customary to define also the 
hyperuniversal combination 

Rj ^ (A-y/'f^, (A19) 

where is the amplitude of a transverse correlation length defined from the stiffness 
constant p^, i.e. = Ps^- Rj can be determined directly from experiments below (see 
e.g. Ref. [0). 



APPENDIX B: GENERAL DISCUSSION ON THE PARAMETRIC 

REPRESENTATIONS 



A wide family of parametric representations was introduced a long time ago p!^-P 
in forms that can be related to our Eqs. (|33|) . Since the application of parametric repre- 
sentations in practice requires some approximation scheme, one may explore the freedom 
left in these representations and understand how this freedom may be exploited in order to 
optimize the approximation. The parametric form of the equation of state forces relations 
between the two functions m = pm{6) and h{6), but it is easy to get convinced that one of 
the two functions can be chosen arbitrarily. For definiteness, let's take m to be arbitrary 
and find the constraints that must be satisfied by h{6) as a consequence of the equation of 
state. 

It is convenient for our purposes to establish these constraints by imposing the formal 
independence of the function F{z) from the parametrization adopted for m = pm{6), which 
we may symbolically write in the form of a functional equation, 



_5_ 

5m 



ph{e) 



0, 



(Bl) 
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keeping z fixed. By expanding m[ff) according to Eq. (|38D and treating the coefficients p 
and c„ = m2„+i as variational parameters, we may turn the above equation into a set of 
partial differential equations (keeping z fixed) 



d_ 

dp 
d 

dci 



ph{0) 

(1 - e2)/3<5 



0, 
= 0, 



(B2) 
(B3) 



which must be satisfied exactly for all i by the function h{9). Simple manipulations lead to 
the following explicit form: 



Y{9) 



m{9) 



(B4) 



Bh 

OCi 



l-9')^ + 2(3S9h 



(B5) 



where Y{9) is defined in Eq. (^). In turn, by expanding 

oo 

hi9,p,c,)=9 + J2hn+i{p,c,)9^''^\ 



(B6) 



n=l 



and substituting into the above equations one obtains an infinite set of linear differential 
recursive equations for the coefficients /i2n+i; which generalize the relations found in Ref. []TD|, 
where the case m{9) = 9 was analyzed. A typical approximation to the exact parametric 
equation of state amounts to a truncation of /i to a polynomial form. We may in this 
case refine the approximation by reinterpreting the ffist recursion equations involving a 
coefficient /i2i+i which is forcefully set equal to zero as stationarity conditions, which force 
the parameters p and Cj into the values minimizing the unwanted dependence of the truncated 
F{z) on the parameters themselves, i.e. on the choice of the function m = pm{9). The above 
procedure implies global (i.e. 6'-independent) stationarity, and as a consequence all physical 
amplitudes turn out to be stationary with respect to variations of p and q. 

These statements are fairly general, but it is certainly interesting to consider the ffist few 
non trivial examples. For the lowest order truncation of h we may adopt the parametrization 




hie) 



which includes both the Ising model in three dimensions [p 
models with Goldstone bosons in d dimensions [p = 2/{d ^ 
the following relationship must then hold: 



6 



p + ci = 7 - 



p_ 
01 



(B7) 

1) and general 0{N) symmetric 

2) ). It is easy to recognize that 



(B8) 
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Global stationarity implies that the stability conditions may be extracted from the variation 
of any physical quantity. In particular we may concentrate on the universal zero of F{z), 
Zq, noting that Zq = -2(6*0), cf. Eq. (|36|). As a consequence of the above results, the simplest 



models can all be described by the parametrization 



zo_ ^ ^{i-ci)ei-p 
V6 (1 - eiY ^ 



ciOt). (B9) 



Let us first consider the case ci = 0. The minimization procedure leads to 

2 ^ 67(7 -p) 
P 7 - 2p/5 ' 

Setting p = 1 one immediately recognizes the linear parametric model representation for the 
Ising model |[T^n|. Unfortunately when p = 2 the solution is not physically satisfactory. 



because it gives 6^ < for all > 2, and therefore the scheme is useless for models with 
Goldstone singularities. 

Let us now include ci. Requiring zo to be stationary with respect to variations of both 
parameters, we obtain 



P 



2 



67(7 - P + 1) 
37-2(p-l)/5' 
37 - 2{p - 1)(3 

(3-2/5)7 ' 
2(7 -p) + (2/3-1) 

37-2(p-l)/5 ' ^^^^^ 



implying also 



- /3 / . \P 



(B12) 



V6 \ 3-2(3 J \2(3^ 
In the Ising model the above solution reduces to 

= 27, (B13) 



° 3-2/3' 
ci = ^(/5 + 7)-l. 

Note that, substituting the physical values of the critical exponents (3 and 7 for the = 1 
model, Ci turns out to be a very small number (ci = 0.04256) and the predicted numerical 
value of Zq is 2.8475, consistent within 1% with the linear parametric model prediction [llO[ . 



It is fair to say that in the Ising case the above solution has a status which is comparable 
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to the linear parametric model, both conceptually and in terms of predicting power. It is 
therefore possible to take it as the starting point of an alternative approximation scheme 
whose higher-order truncations might prove quite effective. 

Unfortunately when we consider the XY system, setting p — 2 and choosing the values of 
the exponents pertaining to N — 2, the value of ci becomes too large for the approximation to 
be sensible. Indeed we get ci = —0.6762 and all testable predictions turn out to be far away 
from the corresponding physical values. It is however worth exploring the features of this 
approach because, as we shall show, it has formal properties which might prove useful when 
considering parametric representations of the equation of state for higher values of N. Let 
us indeed consider the function g2{0) entering the parametric representation of the magnetic 
susceptibility. We know that this function will in general show singularities in the complex 
9 plane corresponding to the zeroes of the function Y{9). However, when substituting the 
expressions of h{6) and m{9) obtained from the saddle-point evaluation of the parameters 
p, 9q and Ci, after some simple manipulations, we find out that all singularities cancel and 



This fact was already observed in the case ci = for all values of the truncation order t. 
Therefore, the stationarity prescription is a way to ensure a higher degree of regularity in 
the parametric representation of thermodynamic functions. 

Finally, let us observe that the stationary solution can be applied to the large- limit 
of 0{N) models in any dimension 2 < o? < 4. In this limit /?=|,7=p = 2/((i — 2). As 
a consequence we obtain from our previous results — 12/ [d -\-2), 9Q — {d-\- 2)/4, ci = 0, 
implying also 



We therefore obtained an exact parametrization of the equation of state in the large- iV limit 
for all d. Thus, for sufficiently large values of N, the scheme we have defined may be a 
sensible starting point for the parametric representation of the thermodynamical functions 
in the critical domain. 







4-d 



(B15) 
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